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Abstract 

We perform a model-independent study for top-antitop and top-top resonances in 
the dilepton channel at the Large Hadtron Collider. In this channel, we can solve the 
kinematic system to obtain the momenta of all particles including the two neutrinos, and 
hence the resonance mass and spin. For discovering top-antitop resonances, the dilepton 
channel is competitive to the semileptonic channel because of the good resolution of lepton 
momentum measurement and small standard model backgrounds. Moreover, the charges 
of the two leptons can be identified, which makes the dilepton channel advantageous for 
discovering top-top resonances and for distinguishing resonance spins. We discuss and 
provide resolutions for difficulties associated with heavy resonances and highly boosted 
top quarks. 



1 Introduction 



Being the unique standard model (SM) fermion with a mass of the electroweak symmetry 
breaking scale, the top quark may be closely related to the TeV scale new physics. In particular, 
many of the new physics candidates predict a tt (tt) resonance, i.e., a heavy particle that 
decays to tt {tt). The tt resonance occurs, for example, in Technicolor [Ij, Topcolor [2j, 
Little Higgs ^, and Randall- Sundrum (RS) models [1], while the tt resonance exists in the 
grand unified theory in the warped extra-dimension [5]. Therefore, it is crucial to study ti 
{t t) invariant mass distributions and look for possible resonances at the ongoing Large Hadron 
Collider (LHC), which may provide us the opportunity for revealing the new physics beyond 
the SM . 

The top quark almost only decays to a 6 quark and a W boson. Depending on how the W 
boson decays, events with a pair of tops can be divided to the all-hadronic, the semileptonic 
and the dilepton channels. The all-hadronic channel, in which both W^s decay hadronically, 
has the largest branching ratio of 36/81, but suffers from the largest background since that all 
observed objects are jets. The semileptonic channel, in which one W decays hadronically and 
the other one decays leptonically, has a significant branching ratio of 24/81 and also smaller 
background. Although there is one neutrino in the event, only its longitudinal momentum is 
unknown, which can be easily extracted using the W mass constraint. Therefore, this has been 
thought to be the best channel for discovering tt resonance and most of existing studies have 
been concentrating on this channel [6l [TJ [8]. The dilepton channel, in which both W^s decay 
leptonically, has been thought to be a very challenging and not promising channel. The reason 
is twofold: first, not counting r's, the branching ratio for this channel is only 4/81; second, due 
to the fact that there are two neutrinos in the final states, the event reconstruction is much 
more difficult than the semileptonic channel. 

Nevertheless, the dilepton channel also has its own merits, making it more than a com- 
plementary to the other two channels. An obvious advantage is that it has much smaller SM 
backgrounds. More importantly, the two leptons in the decay products carry information that 
is unavailable in the other channels. First, it is well-known that the charged lepton is the most 
powerful analyzer of the top spin [9], [10], because its angular distribution is 100% correlated 
with the top polarization in the top rest frame. The down-type quark from hadronic decay 
of the W boson has an equal power, but it is indistinguishable from the up-type quark in a 
collider detector. If the b jet from the top decay is not tagged, the ambiguity is even worse. 
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Only the dilepton channel is free from this ambiguity. 

Secondly, the charges of the two leptons are both measurable, which makes the same-sign 
dilepton channel ideal for studying tt or it production, since it has very small SM backgrounds. 
Note that although we are discussing resonances, the analysis applies equally for any events 
with two same-sign top quarks, as long as there are not missing particles other than the two 
neutrinos. For example, it can be used to study the excess of or tt production in flavor 
violating processes [TTl [T2| IT3] . On the contrary, the charge information in the other two 
channels is unavailable!, and hence a more significant event rate is needed to see an excess over 
the SM t i background. 

Motivated by the above observations, we perform a model independent study on tt (tt) 
resonances in the dilepton channel. The crucial step of this analysis is the event reconstruction, 
which we describe in the next section. We will focus on the most challenging case when the 
resonance is heavy (> 2 TeV) and discuss a few related difficulties and their solutions. As an 
illustration, the method is applied to a KK gluon in the RS model with a mass of 3 TeV. In 
Section [3l we estimate the discovery limits of representative resonances with different spins. 
It is shown that despite the smaller branching ratio, the discovery limits from this channel 
compete with those from the semileptonic channel. In Section HI we present the method for 
spin measurements and estimate the minimal number of events needed to distinguish the spin 
of the resonance. Section [5] contains a few discussions and the conclusion. 



2 Event Reconstruction 
2.1 The Method 

In this section, we discuss the method for reconstructing the 1 1 system in the dilepton channel 
at the LHC The process we consider is pp ^ U ^ ti ^ hbW^W" —* hM^l^v^vp^ with 11 a tt 
resonance and i = e , fi. There can be other particles associated with the 11 production such 
as the initial state radiation, but in our analysis it is crucial that the missing momentum is 
only from the two neutrinos. The method described in this section can also be applied to tt 
resonances. 

Assuming tops and W^s are on-shell and their masses are known, the 4-momenta of the 
neutrinos can be solved from the mass shell and the measured missing transverse momentum 

^It is possible to identify the charges of the &-jets but only at a few percent level. 
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pI = pI = 0, 

{Pv+Pi+f = {Pv+Pi-f = rn^, 
{Pv+Pl++Pbf = {Pu+Pi-+Pbf = 

Pl+Pt=f. Pl+Pl=f, (1) 

where Pi is the four-momentum of the particle i. We have 8 unknowns from the two neutrinos' 
four-momenta and 8 equations. Therefore, Eqs. ([T]) can be solved for discrete solutions. This 
system can be reduced to two quadratic equations plus 6 linear equations [151 [16]. general, the 
system has 4 complex solutions, which introduces an ambiguity when more than one solutions 
are real and physical. After solving for p^, and pp, it is straightforward to obtain pt and pt and 
calculate the ti invariant mass = {pt +Pi)'^- 

The system in Eqs. ([1]) has been applied to measure the top mass [HI ITT] and to study the 
spin correlations in decays [T7l[T8]. These studies focus on low center of mass energies below 
1 TeV and involve only the SM tt production. We will concentrate on the heavy-resonance 
case when t, i and their decay products are highly boosted. There are a few complications in 
disentangling new physics contributions from the SM, as discussed below. 

The first complication comes from the fact that for a highly boosted top, its decay products 
are coUimated and therefore are difficult to be identified as isolated objects. In other words, 
all decay products of the top, in either the hardronic or the leptonic decay channels, form a 
fat "top jet". This interesting fact has triggered recent studies for developing new methods to 
distinguish top jets from ordinary QCD jets [ISl[2nj- For the dilepton channel, in order to keep 
as many signal events as possible, we include both isolated leptons and non-isolated muons. 
Non-isolated muons can be measured in the muon chamber, while non-isolated electrons are 
difficult to be distinguished from the rest of the jet and therefore not included in our analysis. 
This is very different from the low center-of-mass energy case where two isolated leptons can 
often be identified. 

Once non-isolated muons are included, we have to consider the SM non-tt backgrounds 
such as 6 1 and b b productions with one or two muons coming from 6 or c hadron decays. Since 
muons from hadronic decays are relatively softer, we will use a high px > 100 GeV cut for the 
non-isolated muons to reduce the background. This is similar to using the jet energy fraction 
carried by the muon as a cut |J^. Similarly, it is unnecessary to require one or two 6-jet 
taggings, which may have a small efficiency at high energies [21] . Instead, we consider all signal 



constraints [m] : 
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and background events with two high-pT jets. Besides high-p^ cuts, the mass-shell constraints 
in Eqs. ([T]) are also efficient for reducing the background/signal ratio. 

The second complication is caused by wrong but physical solutions. Part of the wrong 
solutions come from wrong combinatorics-either one or more irrelevant jets or leptons from 
sources other than tt are included in the reconstruction equations, or the relevant jets and 
leptons are identified but combined in a wrong way. Even when we have identified the correct 
objects and combinatorics, there can be wrong solutions due to the non-linear nature of the 
equation system. As mentioned before, there could be up to three wrong solutions in addition 
to the correct one. The wrong solutions will change the ti invariant mass distribution. This 
is not a severe problem for a light (< 1 TeV) resonance because both signals and backgrounds 
can be large. The wrong solutions will smear but not destroy the signal peak. For heavy 
resonances in the multi-TeV range, the signal cross section is necessarily small due to the rapid 
decreasing of the parton distribution functions (PDF's). This would not be a problem if we 
only obtained the correct solution since the decreasing would happen for both the signals and 
the backgrounds. However, when a wrong solution is present, it will shift the tt invariant mass 
to a different value from the correct one, either lower or higher. Due to the large cross section 
of the SM 1 1 production in the low invariant mass region, even if a small fraction of masses are 
shifted to the higher region, the signal will be swamped. 

Wrong solutions exist because the momenta of the neutrinos are unknown except the sums of 
their transverse momenta. Clearly, for a 1 1 invariant mass shifted to be higher than the correct 
value, the solved neutrino momenta are larger than their right values statistically. Therefore, 
we can reduce the fraction of wrong solutions by cutting off solutions with unnaturally large 
neutrino momenta. This is achieved by two different cuts. First, we can cut off "soft" events 
before reconstruction. That is, we apply a cut on the cluster transverse mass m^;, defined from 
the measured momenta [7]: 

m^^^ = (^^pl{l+l-bb) + m^{l+l-bb)+^T^ - (^pT{^+l-bb)+f^y , (2) 

where pT{l^l^bb) and m?{l^l^hb) are the transverse momentum and the invariant mass of the 
l^l~bb system, and = \^t\- Second, after reconstruction, we define a cut on the fraction of 
the transverse momentum carried by the neutrinoes, 

r„, = ?i±4<2. (3) 

Pt + Pt 
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As we will see in Section |231 the r^b cut is useful for increasing signal/background ratio. The 
value in Eq. ([3]) is approximately optimized for the examples we consider and taken to be fixed 
in the rest of the article. On the other hand, we choose to explicitly vary the mx^^ cut to 
optimize the discovery significance because it is what the significance is most sensitive to. In 
practice, one could as well optimize all other cuts and obtain better results. 

The third issue is with regard to the experimental resolutions. The smearing of the measured 
momenta modifies the coefficients in Eqs. ([1]). When the modification is small, the correct 
solutions of the neutrino momenta are shifted, but we st,„ obta... real solution J However, 
when the modification is large, it is possible to render the solutions to be complex. Again, this 
effect is more significant when the top is more energetic. The absolute smearings are larger 
(although the fractional resolution is better), which make it harder to have real solutions . For 
comparison, 38% signal events from a 1 TeV resonance have real solutions. The percentage 
decreases to 26% for a 3 TeV resonance. This is based on a semi-realistic analysis detailed in 
the next subsection. 

The best treatment of this problem is perhaps to find the real solutions by varying the visible 
momenta, and then weight the solutions according to the experimental errors. In this article, 
we adopt a much simpler solution, namely, we keep those solutions with a small imaginary part. 
More precisely, we first solve Eqs. ([1]) for p^, and pp. Then we keep all four complex solutions 
and add them to the corresponding lepton and 6-jet momenta to obtain pt and pt- We demand 

\lm{Et)\ < 0.4 \Re{Et)\ , |Im(Ef)| < 0.4 \Re{Et)\ . (4) 

where Et and Ef are respectively the energies of t and i. Similar to the r^,;, cut, the values we 
choose in Eq. (jll) are approximately optimized and taken to be fixed through the rest of the 
article. For events passing the above cuts, we make the 4-momenta of t and i real by taking 
the norm of each component, but keep the sign of the original real part. Note that complex 
solutions always appear in pairs, giving the same real solution after taking the norm. We only 
count it once. 

2.2 Event Generation 

The hard process of — ^ 11 — — > bbi^i'UiUi is simulated with TopBSM [22] in Mad- 
Graph/ MadEvent [25], where 11 denotes the tt resonance. In this article, we will consider a 

^Note that the finite widths of the top quark and the W boson have similar eff'ect, ahhough their 1 — 2 GeV 
widths are neghgible compared with the detector resolutions. 
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spin-0 color-singlet scalar, a spin-0 color-singlet pseudo-scalar, a spin-1 color octet and a spin-2 
color- singlet. The major SM background processes, including ti, bb, cc, bbiu and jj£i, are 
also simulated with MadGraph/MadEvent using CTEQ6L1 PDF's |26]. We choose the renor- 
malization and factorization scales as the square root of the quadratic sum of the maximum 
mass among final state particles, and pts of jets and massless visible particles, as described 
in MadGraph/MadEvent. Showering and hadronization are added to the events by Pythia 
6.4 [27]. Finally, the events are processed with the detector simulation package, PGS4 [28] . 
We have not included theoretical uncertainties in the cross-section calculations, which mainly 
comes from PDF uncertainties at high invariant mass [22]. In Ref. [22] (Fig. 3), it is estimated 
using the CTEQ6 PDF set that the SM ti cross-section has a theoretical uncertainty around 
20%~ 30% at 2 TeV, increasing to about 80% at 4 TeV, which may significantly affect some 
of the results in our analysis. Nevertheless, we note that the PDFs can be improved with the 
Tevatron data [29] at large x, and our focus here is event reconstruction. Therefore, we ignore 
systematic errors in the following discussions. 

The cuts used to reduce the background/signal ratio are summarized below, some of which 
have been discussed in the previous section: 

1. Before reconstruction 

• At least two leptons satisfying: px > 20 GeV for isolated leptons or pt > 100 GeV 
for non-isolated muons. The two highest pt leptons are taken to be the leptons in 
Eqs. (HD; 

• rriii > 100 GeV where mu is the invariant mass of the two highest pt leptons. 

• At least two jets satisfying: px > 50 GeV for 6-tagged, px > 150 GeV for not-b- 
tagged. The two highest px jets are taken to be the b jets in Eqs. ([1]); 

• I^T > 50 GeV; 

• Varying mT^^ cut. 

2. After reconstruction 

. \lm{Et)\ < 0.4 \Re{Et)\ , |Im(Ef)| < 0.4 \Re{Ei)\ ; 

• r,,fe < 2 . 
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The complex solutions are made real using the method discussed in the previous section. There 
can be 0-4 solutions after the above cuts. We discard events with zero solution. For a solvable 
event with n> 1 solutions, we weight the solutions by 1/n. 

2.3 KK gluon as an example 

We illustrate the efficiency of the reconstruction procedure by considering the KK gluon in the 
basic RS model with fermions propagating in the bulk. The KK gluon is denoted by Ilj, which 
has the following couplings to the SM quarks, 

glji = 0.2 Qs, gl = gl = g,, ^^^ = 4^,, g]^ = -0.2g,, (5) 

where gs is the strong coupling constant and q represents quarks in the light two generations. 
With this set of couplings, the KK gluon has a width Fni = 0.153 Mqi, and the branching ratio 
Br{Jl\ ti) = 92.6%. For a KK gluon of mass 3 TeV, the total leading-order cross section 
in the dilepton channel is approximately 10 fb. The parton level m^f distribution is shown in 
Fig. [H together with the SM ti background, also in the dilepton channel. The interference 
between the KK gluon and the SM is small and ignored in Fig. [H Within the mass window 
(Mg - Fg, Mg + Tg) ~ (2500, 3500) GeV, the total number of events is around 770 for the 
signal and 610 for the background, for 100 fb~^. 

Although the SM ti production in the dilepton channel comprises the largest background, 
we have to consider backgrounds from other sources since we are utilizing not-b-tagged jets and 
non-isolated leptons, which can come from heavy flavor hadron decays. The major additional 
backgrounds are: 

1. ti processes in other decay channels, which include the semi-leptonic channel, the all- 
hadronic channel and channels involving r's; 

2. Heavy flavor di-jets, including bb and cc with bb dominating. 

3. Other processes that contain one or more isolated leptons including jjii, bbiv productions. 

The above backgrounds are included in our particle level analysis. In Table [T], we show the 
number of events of the signal and backgrounds before and after the reconstruction procedure. 
The cuts discussed in the previous subsection are applied, with a moderate ttit^^ > 1500 GeV 
cut. Note that these numbers are without any mass window cut, while the kinematic cuts in 
the previous subsection have been applied. Also note that the number of signal events is much 
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Figure 1: The number of events in the dilepton channel of the ti production through a KK 
gluon at the LHC The mass and width of the KK gluon are chosen to be 3 TeV and 459 GeV, 
respectively. The solid (blue) curve is signal+background and the dashed (black) curve is the 
SM ti dilepton background. 

smaller after detector simulation and applying the kinematic cuts: this is because most of the 
leptons are non-isolated when the ti resonance mass is as high as 3 TeV, and we have only 
included non-isolated muons in the analysis. The probability for both leptons from W decays 
to be muons is only 1/4. This fact, together with the kinematic cuts, drastically reduces the 
number of signal events. This reduction also occurs for the SM t i dilepton events with a high 
center of mass energy. 
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Table 1: Number of signal and background events for 100 fb ^ before and after reconstruction. 

From Table [H we can see the effects of the event reconstruction. Before applying the r,^fe 
cut, the reconstruction efficiencies for the signal events and the SM ti dilepton events are ap- 
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proximately equal and around 50%. The efficiencies for the other backgrounds are substantially 
smaller. Moreover, the cut on the variable r,jh, which is only available after the event recon- 
struction, also favors the signal and the SM ti dilepton events. Therefore, we obtain a larger 
S/B at the cost of slightly decreasing significance S/\fB. In the following, we will define the 
significance after the event reconstruction. 

Of course, the effects of event reconstruction are beyond simple event counting. More 
importantly, we obtain the 4-momenta of the top quarks, which are necessary for determining 
the spin of the ti resonance. We will discuss the spin measurement in Section |H We also obtain 
the mass peak on top of the background after reconstruction, as can be seen from Fig. [21 where 
we show the m^f distributions of both background and signal+background for a few different 
m^^; cuts. For the left plot with m-r^j > 1500 GeV, there is a clear excess of events, although 
the mass peak is not obvious. By comparing with Fig. [H we see that S/B is smaller than the 
parton level distribution in the mass window (2500, 3500) GeV ~ (Afc — Fg, Mg + F^). This 
indicates that wrong solutions from the lower irtti background events have contaminated the 
higher mti distribution. As we increase the mx^i cut, the numbers of both signal events and 
background events decrease, but S/ B is increasing, showing that the contamination is reduced. 
The contamination reduction is also confirmed by tracing back the reconstructed m^f to its 
Monte Carlo origin. For the mx^^ cut of 1500 GeV, the reconstructed background mti in the 
mass window of (2500, 3500) GeV is decomposed as: 44% from the SM tt events with original 
mtt smaller than 2500 GeV; 25% from the SM tt events with original mti larger than 2500 
GeV; the other 21% come from other SM backgrounds. The decomposition becomes (in the 
same order as above) {23%, 43%, 34%} for the mr,, cut of 2000 GeV, and {13%, 60%, 27%} 
for the mx^i cut of 2500 GeV. Nevertheless we cannot choose too high a mx^i cut since it can 
reduce S/\fB. For the KK gluon example, the significance is maximized when the m^^; cut 
is around 2000 GeV. More precisely, in the mass window (2500, 3500) GeV for M^, we have 
S/B = 0.69, 1.3, 1.8 and = 4.9, 6.1, 4.5 for my,, > 1500, 2000, 2500 GeV, respectively 

3 Discovery Limits 

Having discussed our strategy of selecting cuts to optimize the discovery limit, we now consider 
the needed signal cross section for different ti and tt resonances at 5cr confidence level (CL) 
at the LHC. 
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Figure 2: The event distributions of tt invariant masses after reconstruction. The sohd (blue) 
histogram is signal+background and the dashed (black) histogram is background only. From 
left to right, we have the cuts on tut , to be 1500, 2000, 2500 GeV. 



3.1 Statistics 

We are interested in heavy resonances with masses from 2 to 4 TeV. In order to cover most of 
the signal events, we examine reconstructed masses in the range of 1.5 — 5.1 TeV. We ignore 
uncertainties from the overall normalization of the signal and background cross sections, which 
can be determined from the low rriu mass events, where the statistics is much better. To utilize 
the shape differences, we equally divide the mass range to A^bin = 18 bins, which amounts to 
having a 200 GeV bin width. In each bin, we define the number of the background events 
as bi, while the number of the signal events as Sj. When the number of events is small, the 
distribution is Poisson. Following [30], we first calculate the Pearson's statistic 



X 



A^bin 



[rii 



(6) 



where = bi + Si is the measured value and Vi = bi is the expected value. Assuming that the 
goodness-of-fit statistic follows a probability density function, we then calculate the p-value 
for the "background only" hypothesis 



P 



1 



2^b.J2r(iVbi„/2) 



(7) 



where A^bin counts the number of degrees of freedom. For a 5 a discovery, we need to have 
p = 2.85 X lO"'' and therefore ~ 65 for A^bin = 18. 

For a particular resonance, we define a reference model with a known cross section. We then 
vary the rriT^, cut from 1.5 TeV to 3.5 TeV in 100 GeV steps to generate different sets of bi and 
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Si. We find the optimized rriT^, cut that maximizes the x^- After optimizing the ttit^^ cut, we 
multiply the number of signal events by a factor of Z to achieve = 65 or 5 cr discovery. This 
is equivalent to requiring the production cross section of the signal to be Z times the reference 
cross section. 

3.2 Discovery limits 

For ti resonances, we choose a representative set of tt resonances with different spins and 
quantum numbers under SU (3) color gauge group. We label spin-0 color-singlet scalar, spin-0 
color-singlet pseudo-scalar, spin-1 color octet and spin-2 color-singlet particles as 11°, 11°, 11^ 
and respectively. For the spin-0 particles, 11° and 11°, we assume that they only couple to 
top quarks with couplings equal to the top Yukawa coupling in the standard model, and hence 
they are mainly produced through the one-loop gluon fusion process at the LHC Their decay 
widths are around 3/(16 tt) times their masses and calculated automatically in the Madgraph. 
For the spin-one particle, 11^, we still use the KK gluon described in Section 1?^ as the reference 
particle, and use the same couplings defined in Eq. ([5]). The decay width of the KK gluon is 
fixed to 0.153 times its mass. For the spin-two particle, 11^, we choose the first KK graviton in 
the RS model as the reference particle, and choose the model parameter, n/Mpi = 0.1, where 
Mpi is the Planck scale and k is defined in [4]. Its decay width is calculated in Madgraph. 

For 1 1 resonances, we study the spin-one particle, which is suggested to exist at the TeV scale 
in grand unified models in a warped extra-dimension [5j. Under the SM gauge symmetries, the 
X gauge boson is the up part of the gauge bosons with the quantum numbers (3, 2, 5/6). It has 
the electric charge 4/3 and couples to up-type quarks with a form gieabcU^l -^IlI^ '^ii + 
{i is the family index; a, 6, c are color indices; C denotes charge conjugate). In general, the 
gauge couplings gi depend on the fermion localizations in the fifth warped extra-dimension. 
However, we do not specify any details of model buildings including how to suppress the proton 
decay, and only focus on the discovery feasibility at the LHC. For simplicity, we model the tt 
resonance the same way as the KK gluon, but flip the sign of one lepton at the parton level. 
We also choose the reference tt production cross section through X as the cross section of tt 
production through the KK gluon described in Eq. (I5l). We fix the decay width of X to be 
10% of its mass. In our analysis, we use the same set of background events as in the tt case. 
There are two main sources of the SM backgrounds for same-sign dileptons. The lepton charges 
from 6-jets can have either sign. Another source is lepton charge misidentifications. There are 
other intrinsic SM backgrounds from processes like ud ^ W^W^du. However, this is a pure 
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electroweak process and hard to pass the reconstruction cuts. We neglect such processes in our 
analysis. 

In Fig. [3], we show values of the multiplying factor Z for the tt{tt) production cross sections 
to have 5 a discovery at the LHC for 100 fb~^ integrated luminosity. Note that we do not change 
the widths of the resonances according to the models described above when multiplying the 
factor Z. Those models serve as reference points only. In obtaining the discovery limits, we 
have ignored all interferences between the resonances and the SM t i productions. As shown in 
Ref. [22], the interference between a KK gluon or a KK graviton and the SM ti productions 
is negligible. For a spin-0 resonance (scalar or pseudo-scalar), a peak-dip structure in the mti 
distribution is generally visible at the parton level if the resonance is produced through gluon 
fusion similar to the SM Higgs. We do not anticipate that the interference will change our 
results significantly. 
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Figure 3: The multiplying factor Z (shown in the figure is its square root) for the production 
cross sections to have 5o" discovery at the LHC with a 100 fb^^ integrated luminosity, as a 
function of ti{tt) invariant masses. 11°, 11°, 11^ and 11^ are spin-0 color-singlet scalar, spin-0 
color-singlet pseudo-scalar, spin-1 color octet and spin-2 color-singlet ti resonances, respec- 
tively. 

The discovery limits for the KK gluon are 3.2 TeV and 3.7 TeV for 100 fb^^ and 300 fb^^ 
integrated luminosity. For comparison, the discovery limit for the KK gluon in the semileptonic 
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channel is 3.8 TeV for 100 fb"^ and 4.3 TeV for 300 fb"^ given in [8]. There they combine the 
invariant mass and top px distributions. If only the invariant mass distribution were used, the 
discovery limit would be reduced by a few hundred GeV. Therefore, the discovery limit in the 
dilepton channel is competitive to the semileptonic channel. Comparing the black (solid) line 
and the orange (thick dashed) line, we have a better discovery limit for the tt resonance than 
the tt resonance when they have the same production cross section. This is because the SM 
background for tt is much smaller than the background for ti. The X gauge boson can be 
discovered with a mass up to 4.0 TeV and 4.4 TeV, respectively, for 100 fb"^ and 300 fb"^ 

If a tt resonance is discovered, it is important to measure the mass. The peak position of 
the rriu distribution in general does not coincide with the true resonance mass, and also shifts 
according to the my^^ cut applied, as can be seen in Fig. [2l We can eliminate this systematic 
error, as well as minimize the statistical error by using the usual "template" method. That is, 
we can generate the rritt distributions for different input masses, and then compare them with 
the measured distribution to obtain the true mass. A detailed study of mass measurement is 
beyond the scope of this article. 

4 Spin Measurements 

The momenta of all particles are known after event reconstruction, which allows us to determine 
the spins of the 1 1 resonances. We first consider the angular distributions of the top quark in 
the 1 1 resonance rest frame. To minimize the effect of initial state radiation, we use the Collins- 
Soper angle [31] defined as the angle between the top momentum and the axis bisecting the 
angle between the two incoming protons, all in the tt rest frame. In the case that the initial 
state radiation vanishes, this angle becomes the angle between the top momentum and the 
beam direction. Using the lab frame momenta, the Collins-Soper angle is conveniently given 



where mu and px are the invariant mass and the transverse momentum of the 1 1 system and 



One can also consider angular correlations among the decay products of top quarks. As 
mentioned in the introduction, the best analyzer for the top polarization is the charged lepton. 



by 





pf, pf are defined by 




(9) 
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Therefore, we examine the opening angle (j) between the direction in the t rest frame and 
the i~ direction in the i rest frame. The parton level distribution for the opening angle has a 
form 

--^ = kl-DcoscP), (10) 
a a cos (p 2 

where D is a constant depending on the ti polarizations, and hence model details. At particle 
level, the distribution is affected by the experimental resolutions and wrong solutions from 
event reconstruction. In Fig. HI we show the particle level distributions of cos 9 and cos (p for 




Figure 4: Distributions of the Collins-Soper angle 6 (left) and the opening angle (right) at 
particle level for different resonances with a mass of 2 TeV and the SM backgrounds. A mass 
window cut (1600 GeV, 2400 GeV) is applied on all solutions. 

4 different tt resonances: a scalar, a pseudo-scalar, a vector boson that couples to left- and 
right-handed quarks equally, and a KK graviton in the RS model. The cuts described in Sec. 12.21 
are applied with mT^^ > ISOOGeV. A mass window cut of (1600 GeV, 2400 GeV) is also applied 
on the solutions to increase S/B. From the left panel of Fig. HI we see significant suppressions 
in the forward and backward regions of cos^, due to the kinematic cuts. Except that, both the 
scalar and the pseudo-scalar have a fiat distribution in cos 6* and are hard to be distinguished 
from each other. The cos 6 distributions for the vector boson and the graviton show the biggest 
difference with respect to each other. As shown in the right panel of Fig. |U the slope of the 
pseudo-scalar distribution in cos0 has an opposite sign to all others, which can be used to 
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identify a pseudo-scalar resonance. Therefore one has to use both distributions to distinguish 
the four particles. 

Given the distributions, we can estimate how many events are needed to determine the spin 
of a tt resonance. We first reform the question in a more specific way: given experimentally 
observed distributions in cos^ and cos0 generated by a particle of spin Sa, we ask how many 
events are needed to decide, at 95% CL, that they are not from a particle of spin Sb. This is 
being done by comparing the observed distributions with Monte Carlo distributions of different 
spins. If the observed distributions are inconsistent with all but one spin, we claim that we 
have identified the spin. Of course, without real data, the "observed" distribution in this article 
is also from Monte Carlo. We quantify the deviation of two distributions from different spins 
Sa and Sb as 

where A^bin is the total number of bins and is equal to 20 by choosing a 0.1 bin size for both 
cos^ and cos0; Ug^^i and rig^^i are the number of events in the z'th bin, which satisfy Yl'^sa,i = 
When = 33, we claim that we can distinguish the spin Sa particle from the spin 
Sb particle at 95% CL, corresponding to the p-value of 2.5 x 10~^, for 19 degrees of freedom 
(here we keep the total number of events fixed, and hence we have one degree of freedom less). 
The number of signal events (after reconstruction) needed to distinguish each pair of spins are 
listed in Table |21 The same cuts as for obtaining Fig. S] are applied. 



Sa\ Sb 


Vector 


Scalar 


Pseudo-scalar 


Graviton 


Vector 




661 (501) 


262 (140) 


316 (122) 


Scalar 


705 (577) 




199 (94) 


771 (455) 


Pseudo- scalar 


275 (182) 


200 (116) 




240 (128) 


Graviton 


356 (243) 


878 (694) 


239 (123) 





Table 2: Number of signal events after reconstruction needed to distinguish a particle of spin 
Sa from spin Sb at 95% CL. The number of background events is fixed to 136, corresponding 
to 100 fb"^ data. All resonance masses are 2 TeV. For reference, the needed numbers of signal 
events without background are given in the parentheses. 

In Table [21 we have shown two sets of numbers. The numbers of events outside the paren- 
theses are the minimum numbers of signal events needed to distinguish the spin for 100 fb~^ 
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data. We use the same cuts as for Fig. HI The number of background events is 136 with the 
ti dilepton events dominating (109). For reference, we also hst in the parentheses the numbers 
of needed events assuming no background. The background distributions are canceled when 
comparing the observed distributions and the Monte Carlo distributions. However, they do 
introduce uncertainties that can significantly increase the number of needed signal events. 

The numbers listed in Table [2] are large but achievable in some models. For example, a 
KK gluon of 2 TeV in the basic RS model yields 230 events for 100 fb~^ in the mass window 
(1600 GeV, 2400 GeV), which is not enough to distinguish it from other spins at 95% CL. With 
300 fb~^ data, we can distinguish it from a pseudo-scalar or a KK-graviton using the dilepton 
channel alone, but will need to combine other channels to distinguish it from a scalar. 

5 Discussions and Conclusions 

An important assumption leading to the fully solvable system is that the only missing trans- 
verse momentum comes from the two neutrinos from the top decays. There are also other 
sources of missing momenta such as neutrinos from heavy flavor hadron decays. But they are 
usually soft and their effects have already been included in the simulation. More challengingly, 
the assumption is invalid when there are other missing particles in the event, for example, 
in supersymmetric theories with R-parity. Consider the process pp tt* ttXiXi the 
minimal supersymmetric standard model, which has the same visible final state particles as 
a ti resonance. We have to be able to distinguish the two cases before claiming a. ti reso- 
nance. Distributions in various kinematic observables are certainly different for the two cases. 
Nevertheless, we find that the most efficient way to separate them is still by using the event 
reconstruction. 

As an example, we have generated 10,000 events in the above MSSM decay chain and let 
both t and i decay leptonically, for = 1500 GeV and m^o = 97 GeV. There are 705 
events which pass the kinematic cuts described in Section [2] with a rriT^^ cut of 1500 GeV. 
Out of those 705 events, only 30 pass the reconstruction procedure, that is, satisfy Eq. (jl]). 
This is to the vast contrast of a resonance of 3 TeV, where a half of the events after cuts 
survive the reconstruction procedure. The difference between the two cases is not difficult to 
understand: for a tt resonance, ignoring initial state radiations, we have t and t back-to-back 
in the transverse plane and their decay products nearly collinear. On the other hand, the 
directions of the two neutralinos in the MSSM case are unrelated and therefore the direction 
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of the missing pt is separated from both of the bi systems. It is then very unhkely to satisfy 
the mass shell constraints simultaneously for both t and i. 

In conclusion, by reconstructing ti and tt events in the dilepton channel, we studied the 
ti and tt resonances at the LHC in a model-independent way. The kinematic system is fully 
solvable from the W boson and top quark mass-shell constraints, as well as the constraints from 
the measured missing transverse momentum. After solving this system for the momenta of the 
two neutrinos, we obtained the t, i momenta and therefore the tt invariant mass distribution. 
The same procedure can also be applied to the 1 1 system. We showed that this method can be 
utilized to discover and measure the spins of ti and tt resonances at the LHC. 

The event reconstruction is the most challenging when the ti resonance is heavy. This is 
not only because of the suppression of parton distribution functions at high energies. More 
importantly, in this case the top quarks are highly boosted and the lepton and the 6-jet from 
the same top decay are highly coUimated. Therefore, the lepton is often not isolated from 
the 6-jet. To solve this problem, we included non-isolated muons, which can be identified in a 
detector. The 6-tagging efficiency may also degrade at high energies, which drove us to consider 
events without 6-tagged jets. In summary, we included all events with two high px (isolated or 
non-isolated) leptons and two high pt (6-tagged or not-tagged) jets passing the kinematic cuts 
described in Section 12. 2[ Accordingly, we have to consider all SM backgrounds containing the 
same final state particles. We simulated and analyzed all major backgrounds and found that 
they can be efficiently reduced by imposing the kinematic cuts and the mass-shell constraints. 

The reconstruction procedure was applied to four ti resonances with different spins. We 
found that despite a smaller branching ratio, the dilepton channel is competitive to the semilep- 
tonic channel in discovering the ti resonance. This is due to the better experimental resolution 
of the lepton momentum measurement and smaller SM backgrounds. For example, the first 
KK gluon in the basic RS model with fermions propagating in the bulk can be discovered at 5 a 
level or better, up to a mass of 3.7 TeV for 300 fb~^ integrated luminosity. We also considered 
the possibility of finding a tt resonance, for which the dilepton channel is the best because it 
is the only channel in which the charges of both tops can be identified. Due to the smallness of 
the SM same-sign dilepton backgrounds, the 1 1 resonance has a better discovery limit than the 
ti resonance. Assuming the same production cross section as the KK gluon, the tt resonance 
can be discovered up to a mass of 4.4 TeV. 

The dilepton channel is also advantageous for identifying the spin of the resonance. We 
considered the top quark angular distribution in the ti rest frame, and the opening angle 
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distribution of the two leptons in their respective top quark rest frames. Combining those two 
distributions, we found that for 100 fb~^ a few hundred signal events (after reconstruction) are 
needed to distinguish the spins of different resonances. 
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